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ABSTRACT 


Future  underwater  missions  will  require  data  transmission  via  satellite.  In 
particular,  the  Office  of  Naval  Research  (ONR)  is  interested  in  experimenting  with 
communications  using  the  GOES  satellite  system,  which  is  government  owned. 
Unfortunately,  communication  antennas  must  point  to  specific  satellites  in  this  system 
and  thus  underwater  vehicles  must  steer  a  specific  course  on  the  surface  during  the 
communication  process.  While  surfaced,  underwater  vehicles  are  subject  to  wind  and 
wave  disturbances  and  it  has  been  suggested  that  control  using  differential  thrust  from 
propellers  may  provide  advantages.  This  thesis  covers  efforts  to  create  and  test  such  a 
steering  autopilot  based  on  the  use  of  the  ARIES  AUV  and  differing  the  voltage  supplied 
to  each  propeller.  It  is  planned  to  use  the  ARIES  in  an  ocean  experiment  to  test  this 
satellite  communication  capability.  This  control  is  embedded  in  the  control  of  ARIES 
during  extended  pop  up  maneuvers  for  GPS  navigational  fixes.  When  surfaced,  not  only 
are  navigational  fixes  obtained,  but  also  data  packets  are  communicated  to  a  command 
center. 
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I. 


INTRODUCTION 


A,  BACKGROUND 

Future  mission  requirements  will  eall  for  satellite  eommunieations  while  the 
ARIES  Autonomous  Underwater  Vehiele  (AUV)  is  surfaeed.  The  Offiee  of  Naval 
Researeh  is  investigating  a  satellite  system  and  eommunieations  paekage  that  would 
support  two-way,  worldwide  eommunieations.  There  are  two  eurrent  satellite  systems 
that  would  be  appropriate  to  these  missions,  the  INTELSAT  and  the  IRIDIUM  Systems. 
Both  systems  are  eurrently  operational  and  would  be  able  to  support  ONR  requirements. 

IRIDIUM  uses  77  satellites  in  Low-Earth-Orbit  (LEO)  to  cover  every  point  on  the 
globe  (Grub,  1991).  The  satellites  are  inexpensive,  lightweight  and  allow  a  data  transfer 
rate  of  2400  baud  (Grub,  1991).  The  IRIDIUM  system  would  support  maritime 
communications  using  L-band  commercial  satellites  and  is  a  good  choice  for  AUV 
satellite  links. 

The  other  system,  INTELSAT,  is  part  of  the  Oceanographic  Data  Link  (ODE)  that 
is  capable  of  worldwide  communications  using  C-band  satellites.  The  link  is  two-way 
capable,  transmits  and  receives  at  900  bits/second,  and  only  requires  an  antenna  the  size 
of  a  quarter  (Gamache,  2001) — ideal  for  a  small,  autonomous  vehicle.  Using  these 
satellites  would  cut  down  on  cost  because  they  require  only  three  satellites  for  full 
coverage  (Gamache,  2001).  ARIES  would  be  able  to  transmit  data  to  a  receiving  station 
many  miles  away  allowing  a  greater  dissemination  of  data  and  operating  range.  Although 
IRIDIUM  has  a  faster  transfer  rate,  ONR  is  interested  in  the  use  of  the  INTELSAT 
system  because  it  makes  use  of  pre-existing,  geosynchronous  government-owned 
satellites. 

Since  the  current  incarnation  of  the  ARIES  Autonomous  Underwater  Vehicle  has 
no  ventral  rudders,  it  has  no  control  while  on  the  surface  save  for  that  provided  by  its 
twin  thrusters.  As  of  the  writing  of  this  thesis,  a  pointing  device  for  positioning  the 
antenna  has  not  been  developed;  hence,  communicating  with  the  satellite  requires  the 
vehicle  to  be  heading  in  a  certain  direction  in  order  to  point  at  the  proper  satellite.  The 
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advantage  of  the  pateh  antenna  is  that  its  broad  beam  width  does  not  require  preeise 
positioning — only  that  ARIES  point  in  the  general  direetion  (Gamaehe,  2001).  The 
vehicle  would  have  to  be  able  to  respond  to  and  compensate  for  such  outside  forces  as 
current  and  wave  action  through  the  use  of  its  thrusters  alone. 


B,  PROBLEM  STATEMENT 

The  work  in  this  thesis  is  aimed  at  determining  whether  a  differential  thrust 
autopilot — responsive  to  waves  and  current — is  able  to  keep  ARIES  within  a  bounded 
course  necessary  for  a  communications  link  to  the  INTEESAT  system  or  whether  a 
ventral  rudder  system  must  be  added.  To  do  this,  experiments  on  the  heading  control  of 
the  surfaced  ARIES  vehicle  were  conducted  using  differential  thrust  between  port  and 
starboard  propulsion  units. 
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II. 


CONTROL  METHODS 


A,  THE  BASICS 

In  order  to  point  at  a  satellite  while  on  the  surfaee  with  the  rudders  completely  out 
of  the  water,  differential  thrust  must  be  used.  A  control  mode  is  necessary  to  regulate  the 
moment  generated  by  the  motors  and  keep  the  vehicle  on  track.  Track  control  is  the 
process  of  keeping  a  vehicle  on  a  designated  track  as  defined  by  a  pair  of  geographic 
waypoints.  Heading  control  is  the  process  of  keeping  the  vehicle’s  course  or  heading  at  a 
prescribed  value.  In  this  work,  the  required  function  is  heading  control.  The  two 
methods  tried  to  maintain  heading  control  were  sliding  mode  control  and  proportional 
control. 

The  control  mode  regulates  the  feedback  of  heading  error  to  determine  an  input  of 
motor  moment.  The  error  follows  the  form 

error  =  current  _  heading  -  track  _  heading  (2.1) 

where  the  error  needs  to  be  bounded  between  ±  n  radians  (180  degrees).  This  distinction 
becomes  critical  as  will  be  seen  later.  Our  goal  when  designing  the  autopilot  was  to  keep 
the  error  within  ±30  degrees,  which  will  still  allow  the  vehicle  to  communicate  with  the 
given  satellite. 

The  track  heading  is  determined  by  taking  the  four  quadrant  inverse  tangent  of  the 
line  drawn  between  two  consecutive  waypoints.  This  heading  is  limited  between  ±180 
degrees  so  that  a  westerly  track  is  actually  denoted  as  -90  degrees.  Since  a  heading  of 
270  degrees  (due  West)  would  always  result  in  a  positive  error  whether  the  vehicle  was  to 
the  right  or  left  of  track,  it  is  of  great  importance  that  the  error  be  limited  to  ±  180  degrees 
and  thus  keep  the  vehicle  turning  in  the  proper  direction  to  correct  its  heading. 

B,  SYSTEM  OF  EQUATIONS 

The  system  of  equations  follows  that  of  a  basic  six  degrees  of  freedom  model  for 
an  underwater  body.  The  x  and  y-axes  are  defined  as  North  and  East,  respectively  while 
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the  z-axis  is  down.  Since  the  focus  is  on  surface  operation,  the  z-axis  is  ignored  along 
with  the  pitch  and  roll  equations  leaving  only  v  (lateral  velocity  or  sway),  r  (turn  rate), 
and  psi  (heading);  the  model  is  taken  from  Healey  (2003): 


1,  Three  State  Model 

mv,  =  -mU^r  +  Y.v^  +  7^  +  Y.r  +  Y^r  +  Ygd^{t) 

+  N,r  +  N/  +  Ngd^  {t) 


Since  the  rudders  are  out  of  the  water,  Y5,  N3,  and  5r(t)  in  (2.2)  are  0.  A  new 
input  term,  5e(t),  is  needed  to  describe  the  differential  moment  from  the  motors.  The 
added  mass  and  force  coefficients  used  in  (2.2)  were  taken  from  Johnson  (2001). 

Equation  (2.2)  in  matrix  form  modified  to  measure  thruster  moment  and  account 
for  wave  action  looks  like  the  following: 
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(2.3) 


This  matrix  follows  the  basic  form 


Mi:  =  Ax  +  Bm +  FJ  (2.4) 

•  M  is  the  mass  matrix 

•  A  is  the  vehicle  dynamics  matrix 

•  B  is  the  input  parameters  matrix 

•  u  is  the  motor  input 

•  F  is  the  disturbance  matrix 

•  d  is  the  wave  disturbance. 


For  the  purpose  of  initial  study,  the  waves  were  modeled  as  a  sine  wave  with  a 
two  second  period,  and  coefficients  were  estimated  to  match  the  experimentally 
determined  course  change  of  twenty  degrees. 
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disturbance  = 


sin(;z-0 


(2.5) 


The  waves  produce  two  effects  on  the  vehicle:  a  lateral  force  (Y^)  and  a  moment  (Nw). 


C.  SLIDING  MODE  CONTROL 

Sliding  mode  control  follows  the  form  (Healey,  2003) 

u  =  -kx  -  Tjsat  sgn(— )  (2 . 6) 

<!> 

where  u  is  the  moment  imparted  upon  the  vehicle  with  units  of  Nm.  Poles  were  chosen  at 
0,  -0.2,  and  -0.4  rad/sec  to  allow  for  a  slow  vehicle  response,  and  the  gain,  k,  was  found 
using  the  MATLAB  ‘place’  command  (refer  to  the  Appendix  for  the  code  used).  The 
closed  loop  form  (Ac=A-Bk)  was  found  from  the  A  and  B  matrices  from  equation  (2,3)). 
The  eigenvalues  of  Ac  determined  s  such  that  A^s  =  0 .  Finally,  the  sliding  surface,  a, 

was  found  from  s^x  where  x  =  x  -  (Healey,  2002).  With  the  calculated  values 

substituted  in.  Equation  (2.6)  takes  the  form: 

u  =  7.8 124v - 8.5766r -\0sat sgn(-0.0842v  +  0.8722r  +  0.4708^i^)  .  (2.7) 

Since  the  vehicle  cannot  control  v  (sway),  the  equation  simplifies  to 

u  =  -8.5766r-105atsgn(0.8722r  +  0.4708^ii^).  (2.8) 

The  saturation  limit,  p,  was  set  to  10  to  limit  the  moment  to  within  what  the  thrusters 
could  reasonably  produce. 

Refer  to  Figure  1  for  the  SIMULINK  drawing  of  the  system  using  sliding  mode 
control.  SIMULINK  was  used  to  model  the  system  and  test  for  initial  stability  before  the 
vehicle  was  taken  out  on  the  water 
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Figure  1  SIMULINK  diagram  of  sliding  mode  system 


D,  PROPORTIONAL  CONTROL 

The  second  method  tried  was  proportional  control.  A  simple  proportional 
feedback  loop  was  built  using  a  hyperbolic  tangent  as  a  switching  term  to  fdter  out  wave 
noise.  The  control  law  was: 

M  = -1.3*tanh(30*^ii^)  (2.9) 

The  maximum  moment  is  ±  1.3  Nm  and  occurs  approximately  6  degrees  off  desired 
heading. 

Figure  2  contains  the  SIMULINK  drawing  of  the  proportional  system.  Wave 
disturbance,  F,  is  modeled,  and  each  state  can  easily  be  monitored  by  use  of  the  various 
scopes.  Again,  SIMULINK  was  used  to  approximate  the  system  stability  and  test  to  see 
that  the  control  was  operating  in  the  correct  direction  to  correct  heading. 


6 


heading  errorl 


Figure  2  SIMULINK  diagram  of  proportional  control 


Two  control  modes,  sliding  mode  and  proportional,  both  making  use  of 
differential  thrust  were  used  to  try  and  steer  ARIES  on  the  surfaee.  The  following 
seetions  doeument  the  experimental  results  and  eompare  the  effeetiveness  of  both  eontrol 
modes. 
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III.  EXPERIMENTAL  TESTING  AND  RESULTS 


A.  PROCEDURE 

The  sliding  mode  and  proportional  eontrol  eodes  diseussed  in  the  last  ehapter 
were  implemented  in  the  ‘exeef.e’  fde,  and  added  to  the  portion  of  the  eode  exeeuted 
when  ARIES  eame  to  the  surfaee  to  get  a  GPS  update — the  ‘GPS  popup’  seetion.  In  the 
aetual  setup,  ARIES  eould  update  its  position  as  it  transmitted  information  via  the 
satellite  link.  Refer  to  the  Appendix  for  the  segment  of  eode  used  to  eontrol  the  thrusters. 

Testing  involved  an  extensive  experimental  setup.  The  AUV  group  consisted  of 
two  sections:  a  land  group  and  a  sea  group.  Eand  setup  included  a  Command  Trailer 
fitted  out  with  radio  communications  and  the  laptop  computers  used  to  program  and 
monitor  the  vehicle.  Although  the  vehicle  was  completely  autonomous,  program  changes 
could  be  made  in  the  command  trailer  and  radioed  via  the  radio  modem  to  ARIES 
without  the  need  for  bringing  the  vehicle  back  to  base  for  reprogramming.  Program 
changes  and  track  updates  could  be  made  in  response  to  any  problems  or  changes  that 
came  up. 

The  maritime  portion  of  the  project  consisted  of  a  Boston  Whaler-class  vessel 
equipped  with  radio  antennae  for  the  modem  and  VHP  communicator.  The  whaler  would 
tow  ARIES  to  the  operating  area  and  follow  her  during  runs.  The  whaler  crew  could 
watch  the  vehicle’s  performance  and  be  on  hand  for  any  emergencies. 

The  missions  followed  a  box  pattern  created  in  the  ‘Track.out’  file  with  the 
vehicle  coming  to  the  surface  for  GPS  readings.  The  ‘Track.out’  file  simply  contained 
the  geographic  waypoints  that  ARIES  would  follow  during  its  mission  along  with  other 
operational  parameters  such  as  length  of  GPS  reading  and  cruising  depth.  During  these 
‘GPS  pop  ups’,  differential  thrust  control  would  be  active  and  the  heading  measured 
against  the  planned  track  in  order  to  determine  the  error.  Each  specific  pattern  is 
described  in  the  following  sections.  All  operations  were  conducted  in  the  Monterey  Bay. 
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B.  RESULTS 

1.  31  March  2003 

The  first  run  used  a  300-meter  North  leg,  a  50-meter  East  leg,  and  300-meter 
South  leg.  GPS  pop  ups  were  seheduled  for  the  North  and  South  legs  and  were  35 
seeonds  long. 


heading 


Figure  3  31  March:  heading  of  entire  run  #1 

Figure  3  shows  a  graph  of  heading  versus  time  for  the  first  run.  The  first  portion 
(from  0  to  50  seconds)  shows  the  turn  as  the  vehicle  heads  to  the  first  waypoint,  and  each 
horizontal  line  represents  the  north  and  south  legs. 
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Figure  4  31  March:  North  run  #1;  showing  number  of  satellites  in  blue  and 
heading  in  green 

Figure  4  shows  a  close  in  view  of  the  north  run.  The  blue  line  indicates  the 
number  of  satellites  seen  by  ARIES  and  is  greater  than  zero  only  when  ARIES  is  on  the 
surface.  Erom  time  200  to  time  210,  ARIES  is  surfaced  and  remains  on  a  course  of -20 
degrees  or  340  degrees  true.  This  is  the  course  necessary  for  ARIES  to  reach  the  first 
waypoint. 
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Figure  5  31  March:  South  run#l;  satellites  (blue)  and  heading  (green) 


Figure  5  shows  the  southerly  portion  of  the  first  run  and  again  shows  a  steady 
course  of  190  degrees. 
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Figure  6  31  March:  run#l;  thrust  effort  between  port  and  starboard  screws 
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The  thrust  effort  is  shown  in  Figure  6;  note  there  is  no  separation  between  port 
and  starboard  thrusters.  The  control  law  should  have  generated  a  thrust  differential  and 
the  fact  that  it  does  not  points  to  an  error  in  the  code. 

2.  9  April  2003 

Changes  were  made  to  the  code  to  amplify  the  thrust  differential.  The  first  run 
followed  the  setup  for  31  March  with  the  exception  that  the  GPS  runs  were  70  and  100 
seconds  long.  Table  1  includes  a  summary  of  all  data  from  runs  made  on  9  April  2003. 


thrust  effort 


Figure  7  9  April:  Run  #1;  red  port  thrust  (Volts),  green  starboard,  blue 
number  of  satellites. 

Figure  7  now  shows  good  separation  between  port  and  starboard  motors  as 
expected. 
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Figure  8  9  April:  run  #1  heading 

In  Figure  7,  the  port  thruster  is  ahead,  and  Figure  8  shows  the  vehicle  turning  to 
starboard.  The  vehicle  is  to  the  right  of  course,  and  the  thrusters  should  generate  a 
negative  moment  (starboard  ahead)  to  correct.  It  turns  out  that  since  the  heading  is 
negative,  and  the  error  is  determined  as  the  difference  between  heading  and  track  (see 
equation  (2.1)),  the  error  will  always  be  negative — no  matter  whether  ARIES  is  to  the 
right  or  left  of  track.  A  negative  error  will  always  generate  a  positive  moment,  and  the 
vehicle  will  always  turn  to  starboard  as  is  seen  in  the  figure.  The  heading  measurement 
was  corrected  in  the  May  run  as  the  problem  was  not  discovered  until  then.  The 
solution  is  to  limit  the  error  between  ±  n  radians  so  that  the  vehicle  always  turns  the 
shortest  distance  to  attain  the  proper  heading. 

a.  Second  Run 

The  second  run  is  set  up  in  the  same  manner  as  the  first  and  is  included  for 

repeatability. 
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thrust  effort 


Figure  9  9  April:  run  #2;  red  port 

In  Figure  9,  there  is  again  good  separation  between  thrusters  showing  that  the 
eontrol  seheme  is  putting  forth  effort  to  eontrol  the  vehiele’s  heading. 

In  the  southerly  portion  of  Figure  10,  the  vehiele  is  turning  to  starboard  with  the 
starboard  thruster  ahead — the  thrust  is  not  enough  to  overeome  the  wave  aetion  pushing 
against  the  vehiele,  and  ARIES  turns  nearly  180  degrees. 

b.  Fourth  Run 

The  fourth  run  eonsisted  of  four  legs  starting  north  for  300  meters  then 
east  for  500  meters,  south  for  50  meters  and  then  baek  west  for  500  meters.  ARIES  eame 
to  the  surfaee  for  30  seeonds  on  the  North  run,  70  seeonds  for  the  east  run,  and  70 
seeonds  for  the  western  traek. 
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Figure  10  9  April:  run  #2  heading 
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Figure  119  April:  run  #4;  red  port 
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Figure  12  9  April:  run  #4  heading 


In  Figure  12,  the  090  leg  has  a  positive  30  degree  eourse  change  on  the  surface 
while  the  western  leg  changes  40  degrees  to  starboard.  Although  Figure  1 1  shows  the 
thrusters  operating  to  correct  heading,  the  thrusters  are  ineffective.  It  is  apparent  that  the 
motors  are  not  strong  enough  to  overcome  wave  action. 

3.  30  April  2003 

The  following  runs  use  the  proportional  controller  developed  in  Chapter  II 
Section  D.  The  vehicle  path  is  a  300  meter  North-South  by  200  meter  East-West 
rectangular  box  with  GPS  pop  ups  of  50  seconds  on  each  leg.  Refer  to  Table  2  for  a 
summary  of  the  data. 
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a.  First  Run 

Although  this  run  was  entirely  underwater  due  to  an  error  in  depth  eontrol,  it 
aetually  brings  to  light  some  useful  data.  When  the  differential  thrust  mode  is  active,  the 
turn  rate  of  the  vehicle  increases  greatly  as  can  be  seen  in  Figure  13.  The  motors 
generate  enough  moment  to  turn  the  vehicle  while  it  is  underwater  and  unaffected  by 
waves. 


turn  rate  and  thrust  effort 


Figure  13  30  April  #1  turn  rate  and  thruster  voltage 
b.  Second  Run 

As  we  were  having  trouble  controlling  depth,  it  was  set  to  zero  for  all  of 
runs  two  and  three.  Because  of  its  depth,  ARIES  picks  up  satellites  for  most  of  the  run 
even  though  there  were  only  four  scheduled  GPS  readings.  The  only  sections  of  interest 
for  this  thesis,  however,  are  the  sections  of  thrust  control  noted  by  the  diverging  red  and 
green  lines  in  Figure  14  and  Figure  18. 
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thrust  effort 


Figure  14  30  April:  run  #2 

Note  for  the  third  GPS  reading  in  Figure  14,  ARIES  did  not  make  it  to  the 
surface  and  thus  sees  no  satellites. 

Figure  15  shows  the  heading  throughout  the  entire  run.  The  east  leg  is 
bounded  within  ten  degrees  of  ordered  course,  which  is  well  within  communication 
parameters,  while  the  west  leg  is  stable  although  offset  by  thirty  degrees.  For  a  closer 
view  of  the  north  run,  refer  to  Figure  16. 
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Figure  15  30  April:  run  #2;  the  heading  is  in  green  while  the  number  of 
satellites  in  blue 
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Figure  16  close  up  view  north  portion  of  second  run 
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As  can  be  seen  in  Figure  16,  the  north  leg  is  very  stable.  Although  it 
seems  very  uneven,  the  small  seale  magnifies  the  wave  action;  the  course  only  varies  4 
degrees  to  starboard  and  4  degrees  to  port  from  the  average  eourse  of  012. 
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Figure  17  30  April:  run  #2 

The  turn  rate  in  Figure  17  is  problematic  as  it  is  almost  always  positive 
even  when  the  starboard  thruster  is  ahead — ^whieh  should  generate  a  negative  turn  rate. 
This  shows  that  the  wave  action  is  too  strong  for  the  thrusters  to  overcome  no  matter 
which  direction  the  vehicle  is  heading. 


c.  Third  Run 

To  ensure  repeatability,  a  third  run  with  the  same  parameters  as  the  second 

was  added. 
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thrust  effort 


Figure  18  30  April:  run  #3 


Again  we  see  in  Figure  18  that  a  large  portion  of  the  run  is  on  the  surfaee, 
but  the  main  points  of  interest  are  the  four  pop  ups  where  differential  thrust  eontrol  takes 
over. 
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Figure  19  30  April:  run  #3 
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On  the  way  to  the  first  waypoint,  ARIES  steers  a  course  of  014  but  varies 
by  16  degrees  as  is  shown  by  the  oscillation  in  Figure  19.  ARIES  maintains  an  average 
course  of  090  and  170  for  the  east  and  south  runs,  respectively,  but  varies  nearly  20 
degrees  for  the  east  run  and  nearly  30  degrees  for  the  south  run.  Heading  is  not 
maintained  for  the  west  run  where  the  maximum  error  is  over  30  degrees  as  is  seen  in 
Table  2. 


heading 


Figure  20  East  portion  of  run  #3 

The  average  heading  in  Figure  20  is  090  degrees,  and  the  course  varies  by 
12  degrees  to  port  and  6  degrees  to  starboard,  which  is  inside  the  acquisition  range  of  the 
C-band  antenna. 
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heading 


Figure  21  South  portion  of  run  #3 

The  southerly  leg  in  Figure  21  averages  170  degrees,  which  is  within 
parameters,  but  it  varies  by  nearly  30  degrees  (10  to  starboard  and  18  to  port)  during  the 
course  of  its  run.  Although  this  is  still  within  the  beamwidth  of  the  antenna,  it  is  a  very 
large  deviation  and  not  indicative  of  good  control. 
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turn  rate  and  thrust  effort 


Figure  22  30  April:  run  #3 

The  turn  rate  is  again  mostly  positive  as  shown  in  Figure  22  despite 
thruster  orders  to  the  contrary.  Once  more,  this  supports  the  conclusion  that  the  thrusters 
are  not  powerful  enough  to  compensate  for  the  wave  action. 

4.  9  May  2003 

For  this  run,  there  were  two  to  three  foot  waves,  and  this  greatly  affected  ARIES’ 
performance.  The  rudders  were  set  to  zero  while  surfaced  in  order  to  isolate  the  action  of 
the  thrusters.  Table  3  contains  the  results. 

In  the  control  code  (‘execf.c’),  the  moment  coefficient  was  changed  to 

delta_newtons=delta_nm*6 . 5617/2.0 

in  order  to  better  reflect  ARIES’  width  of  24  inches  (0.6  m)  The  term  delta jiewtons  is 
the  force  differential  in  Newtons  and  is  defined  as  deltajiewtons  =  portjhrust  - 
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starboard Jhrust.  The  term  deltajim  is  the  moment  (Newton-meters)  generated  by  the 
differential  thrust. 


heading  and  number  of  satellites 


Figure  23  9  May:  heading  for  entire  run  and  number  of  satellites 

Figure  23  shows  the  heading  inereasing  to  starboard  for  every  GPS  popup.  This 
was  due  largely  to  the  ‘wrapping’  error.  Sinee  ARIES  determined  the  heading  to  be 
negative,  the  errors  for  the  north,  east,  and  south  portions  were  all  negative  and  thus  all 
moments  positive.  The  exception  was  the  western  run  where  the  desired  track  was  -90 
degrees.  Here,  the  error  was  positive,  and  the  motors  generated  a  negative  moment.  The 
turn  rate,  however,  was  still  positive  (refer  to  Figure  25)  because  of  the  waves,  and  thus 
the  vehicle  still  turned  to  starboard  though  not  as  much. 
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Figure  24  9  May;  western  track  heading 

Figure  24  shows  the  heading  varying  by  only  12  degrees  during  the  western 
portion  of  the  run,  but  the  course  itself  is  offset  14  degrees  from  the  desired  course  of 
270.  Even  at  maximum  thrust,  the  thrusters  are  not  enough  to  counteract  waves  and 
current;  another  means  of  control  is  required.  The  data  included  in  Table  3  once  again 
shows  that  ARIES  cannot  maintain  heading  while  on  the  surface  in  anything  other  than 
calm  days. 
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Figure  25  9  May:  turn  rate  and  thruster  voltage 

The  first  and  last  GPS  runs  in  Figure  25  show  the  uneontrollability  due  to  waves. 
The  first  seetion  has  the  port  thruster  ahead  and  should  result  in  a  positive  turn  rate  but  is 
instead  negative;  likewise,  the  last  seetion  shows  the  starboard  thruster  ahead  yet  has  a 
positive  turn  rate. 

Analysis  of  the  eontrol  eode  determined  that  due  to  the  eontinuous  wrapping  of 
heading,  errors  arose  when  determining  heading  differential  using  equation  (2.1).  A 
heading  of  -360  (as  in  Figure  23  for  the  north  portion  of  the  run),  for  example,  would 
always  show  a  negative  error  even  if  the  vehicle  was  to  the  right  of  intended  course.  A 
negative  error  would  call  for  a  positive  moment  and  would  drive  the  vehicle  even  farther 
off  course.  This  could  easily  be  fixed  by  ‘unwrapping’  the  heading  and  keeping  it 
between  ±  n  and  is  done  so  in  the  next  run. 
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5.  23  May  2003 

For  the  final  run,  all  rudders  and  planes  were  set  to  zero  during  the  GPS  portion  of 
the  track.  The  depth  was  set  to  3  m  for  entire  run  in  order  to  test  thruster  controllability 
without  any  interference  by  waves. 

It  was  found  that  the  moment  coefficient  was  too  weak  to  steer  the  vehicle  and 
was  changed  back  to  its  original  form 

delta_newtons=delta_nm*6.56 1 7 
This  corresponds  to  doubling  the  moment  gain. 

The  heading  error  discussed  in  the  previous  section  was  ‘unwrapped’  to  ±  n. 
ARIES  would  now  take  the  shortest  direction  in  a  turn  and  would  not  drive  in  a  direction 
opposite  to  where  it  was  supposed  to  go.  Table  4  shows  that  the  results  are  much  more 
manageable. 
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turn  rate  and  thrust  effort 


Figure  26  23  May:  turn  rate  and  thruster  voltage 

This  run  was  mainly  used  to  test  underwater  eontrollability  without  rudders.  With 
no  wave  action  to  interfere,  we  are  trying  to  see  if  the  thrusters  could  steer  a  straight 
course.  If  ARIES  maintained  course,  it  would  prove  that  the  thrusters  and  control  code 
were  working  but  were  not  powerful  enough  to  perform  on  the  surface  and  deal  with 
wave  action.  Figure  26  shows  that  the  turn  rate  increases  dramatically  when  the  thrusters 
come  online  but  hovers  around  0  degrees/sec.  meaning  ARIES  steers  a  straight  course. 
This  is  confirmed  in  Figure  27  as  well. 
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heading  and  number  of  satellites 


Figure  27  23  May:  heading  for  entire  run 

The  heading  varies  slightly  but  is  otherwise  stable  during  eaeh  of  the  differential 
thrust  seetions.  A  close  up  of  each  run  is  examined  to  demonstrate  their  stability. 

Figure  28  shows  the  easterly  heading  stays  around  090  degrees  and  only  varies  10 
degrees  to  starboard  and  6  degrees  to  port — ^well  within  communication  range. 
Unfortunately,  since  this  entire  run  is  underwater  and  is  only  used  to  test  ARIES’ 
responsiveness  to  differential  thruster  thrust  it  does  not  help  satellite  communications. 

The  average  heading  in  Figure  29  is  180  degrees  and  only  varies  by  5  degrees  to 
starboard  and  4  degrees  to  port. 
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Figure  28  23  May:  heading  and  thruster  voltage  for  east  leg;  port  thruster  is 
red,  starboard  is  green 


Figure  29  23  May:  heading  and  voltage  for  south  leg 
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Figure  30  23  May:  heading  and  voltage  for  west  leg 

Figure  28  through  Figure  30  show  that  the  ‘wrapping’  error  has  been  solved.  The 
thrusters  act  properly  to  correct  the  heading — generating  a  negative  moment  when  the 
heading  is  to  the  right  of  track  and  a  positive  moment  when  the  heading  is  to  the  left  of 
track. 

In  this  section,  we  tried  both  sliding  mode  and  proportional  control  modes  to 
achieve  heading  control.  Despite  correcting  the  way  the  error  was  measured,  it  does  not 
change  the  fact  that  the  thrusters  are  too  weak  to  exert  any  control.  ARIES  has  no 
reliable  control  while  on  the  surface  and  subject  to  wave  action.  The  data  conclusively 
proves  that  the  motors  cannot  match  the  force  of  the  waves. 
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C.  TABLES  OF  RESULTS 

The  following  tables  eontain  the  results  of  the  experiments  organized  by  date. 
When  several  runs  were  aeeomplished  on  the  same  day,  they  are  ineluded  in  separate 
eolumns.  Eaeh  run  is  further  broken  down  into  the  four  eardinal  direetions.  The 
maximum  error  of  the  run  was  ealeulated  by  subtraeting  the  smallest  heading  from  the 
largest.  The  standard  deviation  is  found  using  standard  mathematieal  praetiees,  and  the 
average  eourse  is  a  time  average  of  the  vehiele’s  heading.  When  there  is  no  data  for  a 
partieular  heading,  either  that  heading  was  not  part  of  the  experiment  (as  in  earlier  runs 
with  only  North-South  or  East-West  runs)  or  ARIES  did  not  make  it  to  the  surfaee  and 
thus  gathered  no  data. 


Table  1  Data  results  for  9  April  2003 


9  APR 

runl 

run2 

run4 

North 

maximum  error 

55 

14 

in  a  turn 

std  deviation 

16.5 

3.8 

average  eourse 

45.9 

12.5 

East 

maximum  error 

27 

std  deviation 

8.9 

average  eourse 

103.3 

South 

maximum  error 

7 

std  deviation 

2.6 

average  eourse 

202.9 

West 

maximum  error 

37 

std  deviation 

10.8 

average  eourse 

296.7 
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Table  2  Data  results  for  30  April  2003 


30  APR 

run2 

run3 

North 

maximum 

8 

16 

std  deviation 

2.1 

3.9 

average  eourse 

12.1 

14.5 

East 

maximum 

25 

-17 

std  deviation 

5.6 

4.0 

average  eourse 

95.7 

90 

South 

maximum 

submerged 

29 

std  deviation 

6.1 

average  eourse 

169.3 

West 

maximum 

in  a  turn 

33 

std  deviation 

9.1 

average  eourse 

283.1 

Table  3  Data  results  for  9  May  2003 


9  MAY 

run2 

North 

maximum:  37.5 

std  deviation:  9.2 

average:  15.4 

East 

submerged 

South 

maximum:  92.3 

std  deviation:  30.0 

average  eourse:  268.6 

West 

maximum:  12.0 

std  deviation:  3.4 

average  eourse:  289.7 
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Table  4  Data  results  for  23  May  2003 


23  MAY 

North 

maximum:  9.1 

std  deviation:  2.6 

average:  8.2 

East 

maximum:  16.4 

std  deviation:  4.3 

average:  90.9 

South 

maximum:  9.3 

std  deviation:  2.7 

average:  180.8 

West 

maximum:  9.5 

std  deviation:  2.5 

average:  270.3 
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IV.  CONCLUSIONS 


A.  CONTROL  MODE  CONCLUSIONS 

As  seen  in  the  Chapter  III,  sliding  mode  eontrol  was  saturated  by  the  wave  motion 
and  could  not  provide  adequate  control.  In  general,  the  control  scheme  did  not  work  as 
the  vehicle  always  turned  to  starboard  despite  any  control  effort  by  the  motors.  When  the 
thrusters  were  operating  in  the  correct  direction  to  counter  the  vehicle’s  movement,  they 
did  not  provide  enough  control  effort  to  make  any  difference. 

Proportional  control  was  much  better  but  still  not  good  enough  to  overcome  the 
waves.  The  control  code  directed  the  thrusters  in  the  proper  direction  to  counter  the 
waves  and  was  not  saturated  by  the  wave  frequency  due  to  the  fdtering  of  the  hyperbolic 
tangent  function.  The  small  motors,  however,  were  not  enough  to  handle  rough  seas  and 
wave  buffeting. 

It  was  critical  for  both  feedback  methods  that  the  heading  error  be  bounded 
between  ±180  degrees  in  order  to  get  ARIES  to  make  the  shortest  turn.  Without  this 
limit,  ARIES  would  sometimes  turn  to  starboard  even  though  it  was  to  the  right  of  track. 
The  motors  were  acting  to  correct  the  heading,  but  they  were  taking  the  long  way  around 
the  circle. 

When  the  vehicle  was  completely  underwater  and  not  subject  to  wave  action,  the 
differential  control  performed  admirably.  Unfortunately,  satellite  communications  are 
impossible  while  underwater.  The  success  with  controlling  the  heading  using  only 
thrusters,  however,  proved  that  they  could  act  as  a  backup  should  the  rudders  fail. 


B,  RECOMMENDATIONS  FOR  FURTHER  WORK 

It  is  the  assertion  of  this  thesis  that  ventral  rudders  are  necessary  to  control  the 
vehicle  while  on  the  surface.  Another  study  should  be  performed  to  check  the  feasibility 
of  rudders  and  determine  the  control  code  necessary  to  steer  the  vehicle.  An  additional 
method  to  overcome  surface  control  problems  would  be  to  design  and  build  a  movable 
antenna.  This  would  solve  the  problem  of  the  thrusters  being  too  weak  and  place  control 
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responsibility  on  the  motor  turning  the  antenna.  This  motor  would  have  to  eompensate 
for  wave  action  moving  the  vehicle,  but  the  motor  would  not  have  to  be  strong  enough  to 
resist  waves,  only  fast  enough  to  respond  to  them  and  keep  the  antenna  pointing  within 
thirty  degrees  of  the  target  satellite. 

A  further  possible  fix  includes  stronger  motors.  The  motors  currently  mounted  on 
ARIES  can  control  the  vehicle  when  it  is  underwater;  if  they  were  stronger,  perhaps  they 
could  counteract  the  wave  forces  pushing  ARIES  off  course.  An  analysis  of  the  problem 
would  have  to  determine  how  strong  the  thrusters  needed  to  be  and  see  if  it  were  feasible 
to  mount  such  thrusters  on  a  small  underwater  vehicle  such  as  ARIES.  Price  would  have 
to  be  considered  as  well  as  it  might  be  less  expensive  to  install  lower  rudders.  Whatever 
the  method  used  to  control  ARIES  while  on  the  surface,  the  current  method  of  differential 
thrust  is  of  insufficient  strength  to  handle  the  task. 
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APPENDIX 


‘ENGINE5.M’ 

%SI 

rho=1025; 

L=3; 

b=16*.0254; 

T=10*.0254; 

U=1.5;  %  m/s 
area=T*b; 

thrust=0.1*rho*area*  1.5^2;  %N 
width=23*. 0254/2; 
diam=6*.0254;%propeller  diam  in  m 
%Prof  Healey  inputs 
%Jay  Johnson  Model 
%SI 

Yv=-68.16; 

Yr  =  406.3; 

Ydr=0;%  for  thrust  eontrol 
Nv  = -10.89; 

Nr  =  -88.34; 

Ndr=0;  %for  thrust  eontrol 
MY  =  456.76; 

IN  =  215; 

%  initial  wave  estimation 
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Ywave=-.5;  %  Newtons 

Nwave=-.5*(L/2);  %  Newton-meters 

%  wave=sin(pi*t);  %  modeled  as  sin  wave  period  2s 

%mass  matrix 

%  M*xdot=A*x+B*u+F*d 

%xdot=inv(M)  *  A*  x+inv(M)  *  B  *u+inv(M)  *  F  *  d 

%u  is  in  Nm 

M=diag([MY,IN,l]); 

AA=[Yv,Yr,0;Nv,Nr,0;0, 1 ,0]; 

A=mv(M)*AA; 

BB=[0;  1;  0];  B=inv(M)*BB; 

Fl=[l  0;0  1;0  0]; 

F=mv(M)*Fl; 
d=[Ywave;  Nwave]; 

%  [A_dtB_dt]=o2d(A,B,.125); 

%  [A_dt  F_dt]=o2d(A,F,.125) 
eta=l; 

C=[0  0  180/pi];  %oonvert  to  degrees 
Nmax=1.3;  %proportional  eontroller  (Nm) 

%  %slidmg  mode  eontrol 
J=[0  ,  -.2,  -.4];  %rad/s 
K=plaee(A,  B,  J); 

Ael=A-B*K; 

[V,D]=eig(Ael'); 
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s=V(:,3);  %eigenvector  corresponds  to  0  eigenvalue  (D) 

%  u=kl  *(x-xeom)+k2*satsgn(k3(x-xeom)/phi) 

Kl=-inv(s'*B)*s'*A 

%K1(1)=0;  %  no  sway  eontrol 

%K2=-inv(s'*B); 

K2=-10  %adjustable 
K3=s' 

%K3(1)=0;  %  no  sway  eontrol 


‘EXECF.C’ 

(GPS  Portion) 

if(GET_GPS_FIX  -  TRUE) 

psi_errorDIE=psi_oont-SegPsi; 
while(fabs(psi  errorDIF)  >  pi) 

{psi  errorDIF  =  psi  errorDIF  -  dsign(psi_errorDIF)*2.0*pi;} 
deltanm  =  -1.3*dtanh(30.0*psi_errorDIF); 
delta_newtons  =  delta_nm*6.5617; 

if  (fabs(delta_newtons)  >  9.0)  delta  newtons  =  9.0*delta_newtons/fabs(delta_newtons); 
v_rs  =  sqrt((9  -  delta_newtons)/1.0314); 
v_ls  =  sqrt((9  +  delta_newtons)/1.0314); 
if  (v_rs  <  0.0)  v_rs  =  0.0; 
if  (v_ls  <  0.0)  v_ls  =  0.0; 

EeftSorewSpeedControl(v_ls); 

RightSorewSpeedControl(v_rs); 
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